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Abstract 

' The time-dependent Schrodinger equation with time-independent Hamiltonian matrix is a homoge- 

neous linear oscillatory system in canonical form. We investigate whether any classical system that itself 
fH . is linear, homogeneous, oscillatory and conservative is guaranteed to linearly map into a Schrodinger 

^ ' equation. Such oscillatory classical systems can be analyzed into their normal modes, which are mutually 

bJQ, independent, uncoupled simple harmonic oscillators, and the equation of motion of such a system linearly 

' maps into a Schrodinger equation whose Ifamiltonian matrix is diagonal, with h times the individual sim- 

O ' pie harmonic oscillator frequencies as its diagonal entries. Therefore if the coupling-strength matrix of 

* ^ , such an oscillatory system is presented in symmetric, positive-definite form, the Hamiltonian matrix of the 

Schrodinger equation it maps into is h times the square root of that coupling-strength matrix. We obtain 
' a general expression for mapping this type of oscillatory classical equation of motion into a Schrodinger 

^ I equation, and apply it to the real-valued classical Klein-Gordon equation and the source-free Maxwell 

equations, which results in relativistic Hamiltonian operators that are strictly compatible with the cor- 
respondence principle. Once such an oscillatory classical system has been mapped into a Schrodinger 
equation, it is automatically in canonical form, making second quantization of that Schrodinger equation 
QQ ' a technically simple as well as a physically very interpretable way to quantize the original classical system. 

Introduction 
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A time-dependent Schrodinger equation, viewed as a real-valued equation of motion that couples the real 
and imaginary parts of its wave vector, is a homogeneous linear oscillatory canonical classical system (its 
classical Hamiltonian function is the presentation in the appropriate real canonical variables of the quantum 
expectation value of its Hamiltonian matrix). Here we shall see that a general homogeneous linear oscillatory 
conservative classical system's equation of motion can always be linearly mapped into a Schrodinger equa- 
tion, and that this mapping is invertible if the classical system has no zero-frequency normal modes. When 
the oscillatory classical system's coupling- strength matrix is presented in symmetric form and is positive def- 
' inite (i.e., no zero- frequency normal modes), the corresponding Schrodinger equation's quantum Hamiltonian 

matrix comes out to be h times the positive-definite square root of that classical coupling-strength matrix. 
The Schrodinger equation's complex-valued mapped wave vector can be sensibly normalized such that the 
quantum system's energy expectation value equals the oscillatory classical system's energy function. More- 
over, that classical system can then be immediately quantized by means of the very straightforward second 
quantization of the Schrodinger equation that it maps into, which by its nature is in canonical form. This is 
not only technically simple, it is as well automatically accompanied by a detailed physical interpretation — 
e.g., one has a mathematical depiction of classical- wave/quantum-particle complementarity via the linear 
mapping of the original classical oscillatory degrees of freedom (which have Hermitian representation) into 
the second-quantized Schrodinger-equation wave vector's annihilation and creation components (which have 
non- Hermitian representation). Mapping into a Schrodinger equation of the real- valued classical scalar- field 
Klein-Gordon equation with mass parameter m yields a complex-valued scalar wave function and the Hamil- 
tonian operator (|cpp -I- m^c'*)2, which is in accord with the correspondence-principle prescription for a 
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relativistic free particle of mass m [1] . Such mapping of the classical source- free Maxwell equations yields a 
complex- valued transverse- vector wave function and the Hamiltonian operator c|p|, which is relativistically 
appropriate to the massless free photon [2] . 

Oscillatory classical linear systems which are homogeneous and conservative are described by second-order 
equations of motion that have the form, 

q + Kq = 0, (la) 

where K is a, nonvanishing real- valued matrix, all of whose eigenvalues are real and nonnegative, and whose 
real- valued eigenvectors completely span the real- valued vector space on which K naturally operates. Note 
that the use of the terms "vector" and "matrix" in this article is not intended to exclude vectors that have 
a continuum of components (e.g., functions) or matrices that have a continuum of entries (e.g., operators on 
function spaces). However, in the interest of cutting down on notational clutter, all of the didactic generic 
formulas that are presented in this article which involve vector components or matrix entries display only 
the case that these are discrete — that is notwithstanding the fact that the interesting examples which are 
discussed in the last part of this article all have continuum character. 

We note that first- order classical equations of motion which have the simple form, 

s = Ws, (lb) 

also describe homogeneous linear oscillatory conservative classical systems when is a nonvanishing real- 
valued matrix that has exclusively imaginary eigenvalues whose associated complex-valued eigenvectors 
completely span the extended complex-valued vector space on which the real-valued W can operate. This is 
so because Eq. (lb) implies that, 

s - VV^s = 0, (Ic) 

and the matrix —W"^ can be shown to conform to all the requirements stipulated for the matrix K below 
Eq. (la). To see this, note that if s^i is any complex-valued eigenvector of W, with iuj its corresponding 
imaginary eigenvalue, where w is a real number, then because W is real-valued, the particular complex- 
conjugated vector s* is as well an eigenvector of W, but with eigenvalue —iiv. Therefore the real-valued 
vector Su, + s* is an eigenvector of the real- valued matrix —W^ with the real, nonnegative eigenvalue w^. 
In addition, since the complex-valued eigenvectors of W of the form s^j are assumed to completely span 
the extended complex-valued vector space on which W can operate, it is apparent that the real-valued 
eigenvectors of —W^ that have the form s^i + s'^ completely span the real-valued vector space on which 
the real- valued matrix —W"^ naturally operates — and of course the eigenvalues ui'^ of — associated to 
each member of this complete set of its real- valued eigenvectors are themselves real-valued and nonnegative. 
Therefore the nonvanishing real-valued matrix —W^ of Eq. (Ic) possesses all of the properties that are 
required of the real- valued matrix K of Eq. (la). 

It is further to be noted at this point that if the nonvanishing real-valued matrix W is antisymmetric, 
then it automatically fulfills the remaining requirements that are stipulated below Eq. (lb), and, in addition, 
a linear mapping of Eq. (lb) into a Schrodinger equation is immediately manifest. This is so because if W 
is real-valued and antisymmetric, then iW is Hermitian on the extended complex-valued vector space on 
which W can operate. By virtue of its Hermitian property, iW necessarily possesses a complete set of 
complex- valued eigenvectors, for each of which it has a corresponding real eigenvalue. Those real eigenvalues 
of iW correspond, of course, to imaginary eigenvalues of W with the same corresponding eigenvectors, and 
that set of eigenvectors of course completely spans the extended complex- valued vector space on which W 
can operate. In addition, if we multiply both sides of Eq. (lb) by the factor ih, it becomes a Schrodinger 
equation with the Hermitian Hamiltonian matrix ifiW. 

In the next section we shall show that classical equations of motion given by Eq. (la), with the restric- 
tions on the matrix K that are stipulated below Eq. (la), can always be linearly mapped into Schrodinger 
equations — consequently the same is true for classical equations of motion given by Eq. (lb) with the re- 
strictions on the matrix W that are stipulated below Eq. (lb)). This task will be greatly facilitated by the 
fact that the oscillatory classical Eq. (la) can be analyzed into its normal modes, which, of course, behave 
as mutually independent simple harmonic oscillators. It turns out that a classical simple harmonic oscillator 
equation of motion which has the natural angular frequency co can be linearly mapped into a Schrodinger 
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equation for an ultra-basic single-state quantum system whose one-by-one Hamiltonian "matrix" is either 
the real number hu) or the real number —fuv. The classical equation of motion for a collection of such 
mutually independent simple harmonic oscillators (i.e., an oscillatory classical system that has been analyzed 
into its normal modes) correspondingly linearly maps into a Schrodinger equation whose Hamiltonian matrix 
is diagonal, with its diagonal entries corresponding in one-to-one fashion to the angular frequencies of the 
independent simple harmonic oscillators which comprise that particular collection: each such Hamiltonian- 
matrix diagonal entry is a unique one of those angular frequencies times one of the two allowed factors 

±h. 

We now turn to the technical details of the analysis of Eq. (la) into its normal modes, and the subsequent 
linear mapping of such collections of independent simple harmonic oscillators into Schrodinger equations. 

Analysis into normal modes and their mapping into Schrodinger equations 

The real- valued eigenvectors qj of K in Eq. (la) completely span the real- valued vector space on which K 
naturally operates, and each qj corresponds to a nonnegative eigenvalue coj, where we take coj to be real and 
nonnegative. Therefore the qj satisfy eigenvalue equations of the form, 

Kqj = Lu'^Qj (2a) 

It turns out that we can use these eigenvectors qj to construct a matrix S which is invertible, and for 
which the composite matrix S~^KS is in diagonal form, with all of its nondiagonal entries being equal to 
zero, while its diagonal entries embrace all the eigenvalues w| of K. Because of this diagonal form of the 
matrix S~^KS, it will be the case that each of the components {S~^q)j of the transformation S~^q of the 
dynamical vector q of Eq. (la) satisfies an independent simple harmonic oscillator equation whose natural 
angular frequency ojj is the nonnegative square root of one of the eigenvalues uj'j of the matrix K. In short, 
the components of the transformed vector S~^q are the normal modes of Eq. (la). 

We shall now construct the matrix S by filling its columns with the components of a set of linearly 
independent qj , where that set is sufficiently large to completely span the real- valued vector space on which 
K naturally operates, 

S., {q,\. (2b) 

Because the columns of the matrix S are linearly independent and completely span the real-valued vector 
space on which S (and K) naturally operate, the matrix S will have an inverse S^^. In addition, because 
of the eigenvalue equations given by Eq. (2a) and the definition of S given by Eq. (2b), it is readily verified 
that, 

{KS)kj = Ku{qj)i = {Kqj)k = uj]Skj. (2c) 

This result permits us to verify that S~^KS is precisely the diagonal form of the matrix K mentioned below 
Eq. (2a), 

iS-'KS),nj = iS-')mkiKS)kj = io]iS-')mkSkj = Uj]{S-'S)mj = ^]Smj. (2d) 

It is convenient to denote this diagonal form of K as Ks, 

Ks = S-'KS, (3a) 
If we now multiply Eq. (la) through by the matrix and further define, 

qs = S-'q, (3b) 

we obtain from Eq. (la) that, 

qs + Ksqs = 0, (3c) 
which, from Eqs. (3a) and (2d), reads when written in component form, 

dF{qs)j/dt^+UJ-'Aqs)j=Q, (3d) 
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which is a set of mutually independent simple harmonic oscillator equations whose natural angular frequencies 
coj are given by the nonnegative square roots of the nonnegative eigenvalues w| of the matrix K. From Eq. (3d) 
we see that the normal- mode simple harmonic oscillator variables are the components of the vector qs- 

We wish at this point to further linearly map the set of mutually independent simple harmonic oscillator 
equations encompassed by Eq. (3c) into a Schrodinger equation. Notwithstanding that they have the same 
form, Eqs. (3c) and (la) crucially differ in that Ks in Eq. (3c) is known to be diagonal (with real nonneg- 
ative entries on the principal diagonal and uniformly zero entries elsewhere), whereas K in Eq. (la) is not 
guaranteed to be diagonal. When we now attempt to pass to a Schrodinger equation, we obviously do not 
wish to undo the simplicity that having only diagonal matrices present confers on an equation of motion. 
Therefore we now make it a rigid rule; that any attempted further linear mapping of Eq. (3c) into (hopefully) 
a Schrodinger equation may only be attempted with diagonal matrices. This affords an immediate benefit: 
diagonal matrices all mutually commute. 

Now a Schrodinger equation has the form, 

ihtp = H^, (4a) 

and, of course, our cardinal rule stated above requires that the Hermitian matrix H be diagonal. 

Eq. (3c) is second-order in time, whereas the Schrodinger Eq. (4a) is first-order in time. To reconcile this 
difference in order, it is necessary to take V to be a linear mapping of qs, and possibly of qs itself as well. 
Therefore we now make the ansatz, 

^ = iN{Wqs + qs), (4b) 

where the matrices A'' and W are of course both assumed to be diagonal. We make the further assumption 
that the matrix TV is invertible (i.e., has no vanishing entries on its principal diagonal), which implies that 
it simply factors out of the linear, homogeneous Schrodinger Eq. (4a), and therefore is not determined by it. 
From Eq. (3c), we know that qs = —Ksqs- Therefore putting the ansatz of Eq. (4b) into the Schrodinger 
Eq. (4a) results in, 

ihWqs - ihKsqs = HWqs + Hqs, (4c) 

which yields the two equations, 

ihW = H, -ihKs = HW, (4d) 

that have the solutions, 

H = h{Ks)i, W = -iiKs)i, (4e) 

which are consistent with our assumption that H and W arc diagonal matrices, and also imply that H is 
Hermitian. Putting the results of Eq. (4e) into Eq. (4b) together with the definition of Ks given by Eq. (3a) 
and that of qs given by Eq. (3b) yields the desired linear mapping of q and q of Eq. (la) into the Schrodinger 
equation wave vector tjj, and also yields the associated Hamiltonian matrix H of that Schrodinger equation, 

^ = N{{S-'^KS)iS-^q + iS-^q), H = %{S-^KS)^. (4f) 

From Eq. (4f), bearing in mind that both N and {S~^KS)^ are mutually commuting diagonal matrices 
and N is invertible, it can readily be shown that the Schrodinger Eq. (4a) for if) follows from the underlying 
classical Eq. (la) for q. 

We as well note from Eq. (4f) that if all the eigenvalues of K are positive, i.e., the classical system is 
purely oscillatory, then the diagonalized matrix S~^KS is invertible, as is the diagonal matrix {S~^KS)^ , 
and therefore the linear mapping between q and tj: is also invertible. 

An interesting mathematical point is that since the diagonal entries of S~^KS are all real and nonnegative 
(they are the the eigenvalues of K), {S~^KS)i is certainly defined as a diagonal matrix, but multiply so, i.e., 
the signs of the nonzero diagonal entries of {S~^KS)i can be chosen at will. So from a strictly mathematical 
point of view, Eq. (4f) specifies a whole set of distinct linear mappings of the classical q into Schrodinger 
wave vectors ip, with equally distinct Hamiltonian matrices H = h{S~^KS)'^ to accompany each distinct 
linear mapping. 
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Although the Schrodinger Eq. (4a) does not determine the invertible diagonal "normalization" matrix 
A'' of our Schrodinger wave vector tp of Eq. (4f), we can ask if there is an additional physically sensible 
requirement which impinges on the value of that "normalization" diagonal matrix N. 

Now the behavior of quantum expectation values frequently closely parallels that of their classical counter- 
parts, as Ehrenfest's Theorem attests, and that is particularly the case for simple linear systems. Specifically, 
the expectation value of the Hamiltonian matrix H, namely ip*Hip, is a real-valued function of tp and tp* with 
the dimension of energy which is conserved because the time evolution of ip is governed by the Schrodinger 
Eq. (4a) and the Hamiltonian matrix is Hermitian — this conservation of tp*Htp can be explicitly verified. 
The clear classical analog of tp*H'ip is therefore, of course, the classical conserved energy that is associated 
with the Eq. (3c) classical equation of motion. That classical conserved energy is the nonnegative entity, 

£ks (Qs, qs) =' i<isqs + <]sKsqs)/i^l'^), (5a) 

where the dimension and magnitude of the real positive number 7 depends on the dimension and normaliza- 
tion of qs — note that [qsi qs) is required to have the dimension of energy. That £ks [qSj qs) is conserved, 
i.e., that its time derivative vanishes, follows directly from Eq. (3c) itself and the fact that Ks is diagonal. 

Therefore it is completely sensible physically to attempt to determine N by additionally imposing the 
utterly natural requirement that, 

ip*H%p = £Ks{qs,qs) = {qsqs + qsKsqs)/{'2l^), (5b) 

whenever this is possible — we shall see that Eq. (5b) requires the real nonnegative diagonal matrix Kg to be 
positive definite, i.e.. the classical system must be purely oscillatory. Furthermore, the strictly nonnegative 
character of the classical energy Sxsiqs, qs) now precludes the possibility that the diagonal Hamiltonian 
matrix H = U{S~^KS)^ can have anything other than nonnegative entries. Unlike the Schrodinger Eq. (4a), 
Eq. (5b) is, of course, neither linear nor homogeneous in tp. We now reexpress Eq. (4f) in the more compact 
form, 

ij = N{{Ks)hs + iqs), H = n{Ks)i, (5c) 

and substitute the right-hand sides of both the first and second equalities of Eq. (5c) into the left hand side 
of Eq. (5b). For the left and right hand sides of Eq. (5b) to then be able to be equal, the following equation 
involving the diagonal matrices N*, N and {Ks)i must be satisfied, 

N*N{Ks)i =I/{2n-f'^), (5d) 

where I is the identity matrix. Of course this is not possible if {Ks)i has any vanishing or negative diagonal 

entries. If (i^'g) 2 indeed has only positive entries, which can only be the case if the classical system is purely 
oscillatory, then the simplest solution for the diagonal matrix N is one with only real- valued diagonal entries, 
namely, 

N={Ks)-i/{2^^h)^. (5e) 
Putting this determination of N into Eq. (5c) results in the properly normalized Schrodinger wave vector, 

i^ = {{Ks)iqs+i{Ks)-iqs)/i27^h)-2, H = h{Ks)i, (5f) 

which in the more explicit notation used in Eq. (4f) reads, 

^ = {{S-'^KS)iS-'^q + i{S-'^KS)-iS-'^q)/{2^'^h)i, H = h{S-'^ K S)i , (5g) 

where the diagonal matrices S~^KS and {S~^KS)i now both need to be positive definite, and the real pos- 
itive constant 7 comes from the classical energy function £g-ij(g{S~^q, S~^q) of Eq. (5a) that is appropriate 
to the purely oscillatory classical equation of motion system of Eq. (3c), 

£s-^Ks{S-'q, S-^q) = {{S-^q){S-^q) + {S-'q){S-^KS){S-^q))/{2j% (5h) 
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Because the diagonal matrix (S ^KS)i is positive definite, the linear mapping of q and q into tp given in 
Eq. (5g) is invertible, 

q = {{j^n)/2)iS{S-^KS)-i (V' + V*), 9 = -i{{7^n)/2)i S{S-^KS)i {iIj - ^*). (5i) 

While the Eq. (5g) route to the desired invertible linear mapping of Eq. (la) into Schrodinger Eq. (4a) 
is of great generality in principle, in practice it suffers from the need to explicitly know all the eigenvectors 
of K in order to be able to construct S, and, in addition, from the need to explicitly invert S. 

We see from Eq. (5g) that one of the consequences of having the matrix S and its inverse in hand is 
that the Schrodinger equation's Hamiltonian matrix H = ?i{S~^KS)'2 is presented to us in already diagonal 
form. It is certainly not essential that that be the case. In the next section we therefore simply expunge S 
and its inverse from Eq. (5g), which of course will work if K is diagonal. However it quickly becomes clear 
that the resulting expression still works when K is merely symmetric. 

Schrodinger-equation presentation of symmetrically coupled oscillatory systems 

The result of expunging S and from Eq. (5g) is, 

iP = {Kiq + iK-iq)/{2'y^n)i, H = nK^, (6a) 

and if is a real-valued symmetric positive-definite matrix, all the expressions in it still make sense: in those 

circumstances H = HK^ is well defined as a real-valued symmetric positive-definite matrix itself. Therefore 
H is Hermitian, as required, and and are well-defined as real- valued symmetric invertible matrices. 
Furthermore, it is straightforwardly verified that in consequence of the basic oscillatory classical equation of 
motion of Eq. (la), the wave vector ^/j of Eq. (6a) satisfies the Schrodinger Eq. (4a) with the Hamiltonian 
matrix H = HKi given by Eq. (6a). In addition, when K is a. real- valued symmetric positive-definite matrix, 
Eq. (6a) yields, 

rH^ = {qq + qKq)/i2^^), (6b) 
and if 7 has been appropriately selected such that, 

£Kiq,q) = {qq + qKq)/{2^^), (6c) 
has the dimension of energy, then it is clear that, 

LK{q,q) = {qq-qKq)/{2j% (6d) 

also has the dimension of energy. Moreover, it is easily verified that the Eulcr-Lagrange equation which 
follows from the Lagrangian Lpciq, q) of Eq. (6d) is precisely the Eq. (la) classical equation of motion. Now 
the conserved, energy of any classical system that has a Lagrangian L is well-known to be uniquely given by 
{qVqL — L), which, for the particular Eq. (la) case that L is given by Lxiq, q) of Eq. (6d), is straightforwardly 
verified to be Sniq, q), as defined by Eq. (6c). Therefore, Eqs. (6b) and (6c) show that when K is real- valued, 
symmetric and positive definite, then the expectation value of the Hamiltonian matrix which follows from 
Eq. (6a) is equal to the conserved energy of the classical system of Eq. (la), as required. 

Finally, when K is real- valued, symmetric and positive definite, the inverse of the Eq. (6a) linear mapping 
of q and q into tp is readily calculated to be, 

q={{j^n)/2)iK-i{i) + i;*), q = -i{{^^n)/2)iKi{i; - ^j*), (6e) 

which is, as expected, the result of expunging 5* and from Eq. (5i). 

What if the matrix K of the classical Eq. (la) is nonsymmetric? We then first need to find a real- valued 

invertible matrix S such that the similarity-transformed Ks S~^KS is symmetric. Eq. (6a) is extended 

to cover this situation by replacing K by Ks and q by qs '= S~^q, precisely as in Eq. (5f), except that now 
Ks is merely symmetric and positive definite, not necessarily diagonal. 
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In the next section, we use the machinery of Eq. (5f) and its associated Eq. (3c) similarity-transformed 
version of the Eq. (la) classical equation of motion (albeit always bearing in mind that Kg is merely sym- 
metric and positive definite, not diagonal) to show that the real and imaginary parts of xp times the factor 
(2?i) ^ obey a simple first-order coupled equation of motion which can immediately be Hamiltonized and then 
quantized. This second quantization of an oscillatory classical system's linear mapping into a Schrodinger 
equation is a very easy route to that underlying system's quantization, and one which as well automatically 
yields considerable physical insight. 

Hamiltonization and quantization of the Schrodinger-equation presentation 

Taking the real-valued similarity-transformed Ks in both Eqs. (3c) and (5f) to now be, as discussed above, 
merely symmetric and positive definite rather than necessarily diagonal, we note that the real and imaginary 
parts of the wave vector %p of Eq. (5f), each multiplied (for later convenience) by the factor (2?i)i, are given 
by, 

='(?l/2)i(^ + r) = {KsY^qsh, pe = -iih/2)iii;-r) = iKs)-iqs/l, (7a) 
which are readily seen, as a consequence of Eq. (3c), which is a similarity- transformed version of the under- 
lying Eq. (la) classical equation of motion, to satisfy the simple first-order coupled antisymmetrical equation 
of motion, 

qc = {Ks)hc, Pc = -{Ks)hc. (7b) 

With a little effort, it can also be verified that the Eq. (7b) system implies the Schrodinger Eq. (4a) with 

H = h{Ks)i. Moreover, for H = h{Ks)^, where (Ks)^ is real- valued and symmetric, the two equalities 
of the Eq. (7b) system follow from simply the real and imaginary parts of the Schrodinger Eq. (4a). In 
other words, for the situation that we are concerned with here, namely that H = h{Ks)i, where {Ks)^ is 
real-valued, symmetric and positive definite, the real- valued coupled antisymmetrical system of Eq. (7b) is 
completely equivalent to the complex valued Schrodinger Eq. (4a) . 

In addition, Eq. (7b) also follows from a simple bilinear classical Hamiltonian, namely, 

nKsiqcPc) = iqc{Ks)hc + Pc{Ks)hc)/2, (7c) 

via the classical canonical Hamiltonian equations of motion, i.e., 

qc^^p^'HKsiqcPc), Pc = -^qJHKs{lc,Pc), (7d) 

and the fact that {Ks)^ is a real symmetric matrix. 

By putting the definition of {qc,Pc) given in Eq. (7a) into Eq. (7c) we can reexpress our system's classical 
Hamiltonian in terms of its Schrodinger-equation presentation wave vector V' and complex conjugate 

■HKsiqcPc) = irH^ + i>Hr)/2 = rH^, (7e) 

where the last equality in Eq. (7e) follows from the fact that H — h{Ks)i is a real, symmetric matrix. It 
is pleasing to once again see the quantum expectation value of the Hamiltonian matrix H come out to be 
equal to the Schrodinger-equation presentation's classical energy, i.e., to its classical Hamiltonian. 

Since Eqs. (7c) and (7d) assure us that the classical equations of motion of Eq. (7b) obeyed by {qc,Pc) are 
indeed presented in canonical Hamiltonian form, we can now safely quantize this classical system by imposing 
first Dirac's canonical commutation rules on the components of {qc,Pc), and next Heisenberg's equations of 
motion on the now quantized (qcjPc)- Dirac's canonical commutation rules promote the components of 
{QcjPc) into noncommuting Hermitian operators which obey the commutation relations, 

[{Qc)i, iqc)j] = [{Pc)i, {Pc)j] = 0, [{qc)i, {Pc)j] = iTiSij. (8a) 

The Eq. (8a) commutation relations, in turn, imply that the components of the non-Hermitian quantized 
wave vector ip = {qc + ipc) / {2h)i satisfy, in conjunction with the components of this quantized wave vector's 
Hermitian conjugate "0 ^ = (§c — iPc) / (2?i) i , the following commutation relations, 

[ipi , V^,-] = [V^i t , ^ . t] = 0, [ipi , t] = 5y . (8b) 
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These Eq. (8b) commutation relations are the fundamental ones for the components of our Schrodinger- 
equation presentation quantized wave vector, and they give those quantized wave- vector components and their 
Hermitian conjugates, respectively, their well-known interpretation as annihilation and creation operators, 
which is so often crucial to physical understanding. They are as well the key to constructing physically useful 
orthogonal basis sets for the second-quantized Hilbert space that is the result of the imposition of Dirac's 
canonical commutation rules on the components of the dynamical-variable vector {qc,Pc)- 

The Hamiltonian operator for this quantized (i.e., second quantized) Schrodinger-equation presented 
system is obtained by substituting the quantized dynamical- variable vector {qc,Pc) into the system's classical 
Hamiltonian of Eq. (7c), namely by writing down, 

nKsiQcPc) = {qc{Ks)^qc+Pc{Ks)--Pc)/2, (8c) 

which could have ambiguities due to operator- ordering issues, but it is apparent that those do not arise in 
this case. Noting that the quantized dynamical- variable vector {(fcPc) is given in terms of the quantized 
wave vector V' and its Hermitian conjugate by the quantized analog of the two definitions in Eq. (7a), 
namely = (?i/2)^ (^-f ■;/; ^f) andpc = —i{h/2)i{'ip — -tp'^), wc reexpressthe uniquely defined second-quantized 
Hamiltonian operator 'Hks{Qc,Pc) of Eq. (8c) in terms of the quantized wave vector ip and its Hermitian 
conjugate xp^ 

HKA'Qc,Pc) = {ip^Hip + ipHip^)/2, (8d) 

where H = h{Ks)i, a real symmetric positive definite matrix. 

If wc now apply Hcisenbcrg's equation of motion and the commutation rules for the components of the 
quantized -ip and -ip that are given by Eq. (8b) to the second-quantized Hamiltonian operator written in the 
form given by Eq. (8d), we can calculate the time derivative of any component of the Schrodinger-equation 
presentation quantized wave vector ^, 

diPi/dt = {-i/n)[ip,, (ip^Hip + ipHip t)/2] = {-i/h)iiH$)i + {ipH)i)/2 = {-i/n){Hip)i, (8e) 

where the last step reflects the real symmetric character of the Hamiltonian matrix H = h{Ks)i. Thus we 
have shown that, 

ihdip/dt = Hip, (8f) 

i.e., the Schrodingcr Eq. (4a) which the Schrodinger-equation presentation wave vector ip satisfies is also 
satisfied by that wave vector's operator quantization -ip, which itself is, of course, a vector of the annihilation 
operators of the complete set of quantum states which the components of the wave vector xp individually 
describe. 

We next turn to the Schrodinger-equation presentations of specifically the classical Klein-Gordon equation 
and the source-free Maxwell equations. 

The spinless quantum free peirticle from the classical Klein-Gordon equation 

The classical Klein- Gordon equation for the real- valued scalar field <p differs from the classical wave equation 
by a simple mass term [3, 1], 

'(P+i-c^V"^ +oj'^)(P = Q, (9a) 
where to = {{m(P')/%). Eq. (9a) has the form of Eq. (la) with, 

K = -c^W^ +uj'^, (9b) 

which, on the space of real-valued scalar fields, is a real-valued, symmetric, positive-definite operator with 
the dimension of frequency squared. Therefore, starting with Eq. (6a) above and going right through to 
Eq. (8f), we have results that can all be transcribed for the real- valued classical Klein-Gordon equation. 
We need to bear in mind that during this exercise K is specifically defined by Eq. (9b) and that the real- 
valued classical dynamical vector q is defined by the real-valued (p, which, as a real-valued vector, of course 
has a three-dimensional continuous index instead of a discrete one. In such a case the summation that 
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defines index contraction is willy-nilly supplanted by three-dimensional integration, which compels some 
systematic technical changes in the formalism, for example in the dimension of the variables that one deals 
with (summation is over dimensionless indices, integration here involves the three space dimensions) and in 
the fact that Kronecker deltas give way to three-dimensional delta functions. That notwithstanding, most of 
the results properly transcribed to the case of the real-valued classical Klein-Gordon equation remain very 
similar in appearance to the formulas that run from Eq. (6a) through Eq. (8f). 

In particular, Eq. (6a) needs essentially no modification; one simply bears in mind that the operator K is 
given by Eq. (9b), and one replaces the occurrences of q and qhy (j) and (/>. The only remaining issue is one 
of a global reconciliation of dimension, which requires the determination of the parameter 7 that appears 
Eq. (6a) so as to accord with the conventions one intends to adopt for the classical Klein-Gordon theory. Now 
one conventional choice of dimension for ^ is the same as that of the electromagnetic vector potential A [3, 
1], which implies that / |V0|^d^r has the dimension of energy. A glance at the classical conserved energy 
given by Eq. (6c) reveals that 7 must have the dimension of c, so we choose the value c for 7. With that, 
Eq. (6a) yields the mapping into the wave function and Hamiltonian operator of the Schrodinger equation 
that corresponds to the classical Klein-Gordon theory, 

^= {Ki(j, + iK-i^)/{2(^n)i, H = hKi, (9c) 

where the operator K is, of course, given by Eq. (9b). The inverse of this mapping from (j) and (j) into the 
complex- valued Schrodinger wave function is easily calculated, or may be transcribed from Eq. (6e), 

(j) = c{n/2)hK-i{il) + ilj*), ,p=-ic{n/2)iKi{il)-il)*). (9d) 

Now let's take a closer look at the Schrodinger equation's Hamiltonian operator, 

H = hKi =h{-c'V^ + {{mc')/hf)^ (9e) 

In configuration space the quantum momentum operator p is well-known to be given by, 

p = -iW, (9f) 

so that, 

-V^ = |p|V7i^ (9g) 
which, when substituted into the expression for H in Eq. (9e), yields, 

(|cp|2+m2c4)5, (9h) 

which is precisely the quantization of the standard relativistic energy of a free particle of mass m. Thus 
we have the fascinating state of affairs that the classical Klein-Gordon equation (i.e., with real-valued (j)) is 
linearly isomorphic to the very Schrodinger equation with the correspondence-principle mandated square- 
root Hamiltonian for the free particle of mass m that Klein and Gordon were in fact trying to sideline. 
If Klein and Gordon had but been aware of the Eq. (6a) theorem with its square-root character of the 
Hamiltonian matrix H = hK"^ , the history of relativistic quantum mechanics and its second quantization 
might have taken a different route, one in closer harmony with the correspondence principle. 

Second quantization of the Schrodinger wave function ip for the classical Klein-Gordon theory can be 
transcribed from Eqs. (8). Here the different dimension of tp that is imposed by its continuum character 
results in its basic commutation relations coming out in terms of a three-dimensional delta function instead 
of in terms of the Kronecker delta of Eq. (8b). 

[^(r),^t(r')] = 5(3)(r-r'), [V^(r), V^(r')] = 0, [^^ (v),^^ {v')] = Q. (9i) 

This promotion of the Schrodinger wave function '0(r) to operator field is the most straightforward and 
physically transparent route to the quantization of the classical Klein-Gordon field 0(r), which, of course, 
is explicitly given by Eq. (9d) in terms of the Schrodinger wave function and its complex conjugate. The 
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familiar physical interpretation attached to the commutation relations of Eq. (9i) is that the operator field 

■0'l'(r) creates a relativistic spinless particle of mass m at location r, while the operator field V-'(r) destroys such 
a particle. Such particle creation and destruction operator fields are non-Hermitian. Ifowever, from the first 
equality of Eq. (9d) we note that the quantized classical Klein-Gordon field ^(r) itself will, on the contrary, 
turn out to be Hermitian, and will be ambiguously capable of both particle creation and annihilation. In 
light of the second equality in Eq. (9d), the same comments apply to the quantization of the time derivative 
of the classical Klein-Gordon field rf^(r)/rft. A telling characteristic of both of these Hermitian fields is 
that hy themselves they only obey the original second-order real-valued classical Klein-Gordon equation. 
Eqs. (9c), (9d) and (9i) thus mathematically depict the complementarity of the quantized particle outlook 
(oriented toward non-Hermitian second-quantized wave-functions that unambiguously either annihilate or 
create particles, and obey a first-order complex-valued quantum Schrodinger equation) to the classical wave 
outlook (oriented toward Hermitian fields that by themselves only obey a real-valued second-order classical 
wave equation). 

Finally, we wish to exhibit, in terms of these quantized Schrodinger wave functions that create or destroy 
particles, the Hamiltonian operator functional that oversees free relativistic spinless particles in the second 
quantized world (we already met this operator in schematic form in Eq. (8d)), 

#[^,^t] = i/[^t(r)(-c2?i2v2 +TO2c4)5^(r) + V^(r)(-c2;i2y2 ^^2^4)i^t(r)]d3r. (9j) 

We now turn to the similar Schrodinger equation that corresponds to the real- valued homogeneous linear 
source-free Maxwell equations. The difi^erences to the Schrddinger-equation results for the classical Klein- 
Gordon equation are that the resulting relativistic particle is massless, and that its wave function is a vector 
field which is strictly transverse. 

Pree-photon quantum mechanics from the source-free Maxwell equations 

In the source-free case, the Coulomb and Gauss laws tell us that both the electric and magnetic fields are 
purely transverse, i.e., V • E = and V • B = 0. The results of the Maxwell law and Faraday's law in the 
source-free case are, 

E = cV X B, B = -cV X E. (10a) 

This first-order equation system has the simple antisymmetrical character of Eq. (7b), which readily produces 
a Schrodinger equation. For example, the extremely simple transverse-vector wave function ansatz = 
E-|-iB will in consequence of Eq. (10a) satisfy the Schrodinger equation which has the Hamiltonian operator 
he curl. Unfortunately this operator has odd parity, and therefore is not a physically appropriate Hamiltonian 
for electromagnetism. The reason that a Hamiltonian of odd parity has manifested itself here is that the 
transverse vector fields on either side of each of the two equations of Eq. (10a) are of opposite intrinsic parity: 
namely E is a polar vector field, while B is an axial vector field. So it should be feasible to extract a physically 
appropriate even-parity Schrodinger-equation Hamiltonian operator from source-free electromagnetic theory 
by first recasting its linear homogeneous equations of motion such that they involve only transverse vector 
fields which all have the same intrinsic parity. We shall do this here by mapping the transverse axial- vector 
magnetic field B into a transverse polar-vector field that is already well-known to electromagnetic theory, 
namely the vector potential in radiation gauge [4]. Specifically, we define, 

a"^' {-V^)-\VxB), (10b) 

where by (— V^)^^ we mean the standard real-valued symmetric integral operator with the Coulomb kernel. 
Eq. (10b) implies that, 

V • A = 0, (10c) 

i.e., A is a transverse vector field. Furthermore, since B is itself a transverse vector field, Eq. 10b implies 
that, 

V X A = B, (lOd) 
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which is, of course, the basic propcii ty of a vcictor potential A. We can further dehneate the properties of A 
in source- free electromagnetism by using its definition together with Faraday's law (i.e., the second equality 
in Eq. (10a)) to calculate its time derivative, 

A = (-V2)-i(V X B) = -c(-V2)-i(V X (V X E)) = -cE, (lOe) 

where the last equality holds when E is transverse, which is, of course the case for source- free electromag- 
netism. So in that case, 

E = -A/c. (lOf) 

Eqs. (lOd) and (lOf) together imply that for source- free electromagnetism, we can obtain both of B and E from 
A, so we only need to concern ourselves with calculating the polar transverse vector field A. Therefore we 
now substitute Eqs. (lOd) and (lOf) into the Maxwell law, which in the case of source-free electromagnetism 
is the first equality of Eq. (10a), to obtain a linear homogeneous second-order equation which involves the 
polar transverse vector field A alone, 

A-c^V^A = 0. (lOg) 

This is, of coiusc, the classical wave equation, and it bears a m,arked, resemhlance to the classical Klein- 
Gordon equation of Eq. (9a). The only differences are that in Eq. (lOg) the parameter lo that appears in 
Eq. (9a) vanishes identically, and, of course, in Eq. (lOg) the transverse vector field A replaces the scalar field 
(f) of Eq. (9a). Even the d,im,ension of the transverse vector field A is the same as the dimension that we chose 
for (j) by adhering to a common convention [3, 1]. Therefore, for the linear mapping, and its inverse, of the 
real-valued transverse-vector fields A and A into a complex- valued transverse- vector Schrodinger-equation 
wave function we can simply transcribe Eqs. (9b), (9c) and (9d) for the real- valued classical scalar Klein- 
Gordon theory, taking w (and m) to be zero identically, and replacing (j), (j), and by, respectively. A, A, 
and ^. Thus our basic real, symmetric operator is, 

K = -c^V^, (lOh) 

which, to be sure, is not positive-definite in the broadest sense. However, Fourier transformation methodology 
indicates that on a sufficiently restricted function space, — can indeed be regarded as positive definite. 
The operators we actually require in the following mapping formulas are (— V^)^, (— V^)~3 and (— V^)^, 
and they themselves have the tractable-looking positive-definite Fourier representations |k|, |k|~2 and |k| 2 
respectively. 

Transcribing Eq. (9c) as described above, the linear mapping of the real-valued transverse-vector fields 
A and A into the complex- valued transverse- vector Schrodinger-equation wave function ^, together with 
the associated Schrodinger-equation Hamiltonian operator, is given by, 

^ = {KiA + iK-iA)/{2c^h)i, H = nKi. (lOi) 

The inverse of this linear mapping from A and A into the complex-valued Schrodinger-equation wave 
function ^ is, 

A = c(?i/2)*if~4 (* + **), A = -ic{n/2)iKi{^ -^*). (lOj) 

In light of Eq. (lOh) and the fact that in configuration representation p = —ihV, we have from the second 
equality in Eq. (lOi) that the Schrodinger-equation Hamiltonian operator can be written, 

H^hKi = ni-c^V^)i = (c2|p|2)5 ^ c|p1. (lOk) 

This Hamiltonian operator is clearly the quantized version of the relativistic energy of a massless free particle, 
which is appropriate to the free photon, and it as well has even parity. 

By using Eqs. (10b) and (lOe), the vector potential can be removed from the Schrodinger-equation linear 
mapping of Eq. (lOi) in favor of the E and B fields, 

^ = {cK-i{V xB)-iK-iF,)/{2n)i, H = nKh. (11a) 



11 



The mapping of E and B into * given in Eq. (11a) has the inverse, 

B = c(7i/2)5if-i(V X (* + **)), E = i(?i/2)iA'i(* - **). (lib) 

We invite the reader to verify that the complex-valued linear mapping of the classical E and B fields into 
the wave function 'i' which Eq. (11a) specifies, along with its specified Hamiltonian operator H — hKi 
(where K = — c^V^), actually satisfies the Schrodinger equation. (Hint: use the source-free Maxwell and 
Faraday laws of Eq. (10a) and the transverse nature of the source-free E field.) One should also verify that 
the quantum expectation value of the Hamiltonian operator agrees with the classical energy of the E and B 
field system, i.e., that, 

/**(r) • (£r*(r)) (fv = 1/ (|E(r)|2 + |B(r)|2) d^r. (11c) 

Upon their second quantization, Eqs. (11a) and (lib) manifest the expected tantalizing complementary 
interplay of the potential for photon creation and annihilation with the familiar, workaday transverse electric 
and magnetic fields. 

In addition to its zero mass parameter, the second special feature of electromagnetic theory vis-a-vis 
classical Klein-Gordon theory is, of course, the free photon's always transverse polarization (spin) states. 
This signature free-photon characteristic does not cause much in the way of complications, but there is one 
formula concerning second quantization which it notationally impacts, albeit no substantive physical effect 
is involved. The canonical commutation rule for second quantization of the free photon's transverse vector 
wave function might naively be expected to read, 

[(*(r))„(*t(r')),]=<5,,5(3)(r-r'), (12a) 

but this is not mathematically consistent with the transverse character of the second-quantized photon wave- 
functions, i.e., it is mathematically inconsistent with the fact that V • * = 0. The nature of the right-hand 
of Eq. (12a) is one of completeness, but the transverse wave function creation and annihilation operators are 
incomplete in that they do not pertain to vector fields which are the gradients of scalar fields, i.e., they do 
not pertain to vector fields which fail to he transverse. Now the ij components of the projection operator 
onto the subspace of such purely gradient vector fields is given by, 

Pij = -di{-V^)-^dj. (12b) 

We note that Pij is Hermitian, and that its contraction with itself yields itself, which are the two essential 
properties of the ij components of projection operators. Of course its contraction with the components of 
any transverse vector field vanishes. Thus {dij — Pij) are the ij components of the projection operator onto 
the subspace of transverse vector fields, and therefore, 

[(*(r))i,(*t(r')),.] = {rmj-Pij)\r') = (27r)-3 / e*-(-'-') (% -k,k, |k|-2) rf^k. (12c) 

Notwithstanding these fancy maneuvers with projection operators, the only issue which is involved here is 
the simple fact that free-photon creation and annihilation operators (and as well free photon wave functions 
in the first quantized regime) are purely transverse, and therefore a,ny expression involving these operators, 
e.g., the expression which describes their canonical commutation relation, must, of course, correctly reflect 
this fact. There is obviously no physics implication which flows from this requirement of mere notational 
correctness. 

Conclusion 

It is a remarkable fact that any classical system whose equation of motion is linear, homogeneous, purely 
oscillatory and conservative is effectively already first- quantized: once its Eq. (la) coupling-strength matrix 
K has been similarity-transformed to a symmetric, positive- definite presentation, Eq. (6a) invertibly linearly 
maps that equation of motion into explicit time-dependent Schrodinger-equation form with Hamiltonian 
matrix HKi. Thus we see that Michael Faraday and James Clerk Maxwell were actually the first to 
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effectively elucidate a quantized particle, namely the very important and not exactly simple ultra-relativistic 
massless transverse- vector free photon. 

Any complex-valued solution wave function of a time-dependent Schr5dinger-equation has the familiar 
characteristic expansion in terms of the complete set of mutually orthogonal eigenfunctions of that equation's 
Hamiltonian operator. The one-to-one linear mapping of any purely oscillatory linear classical system that 
is homogeneous and conservative into a Schrodinger equation thus implies a characteristic two- component 
eigenfunction expansion of such a classical system's solutions. For the case of certain wave equations that 
fall into the class of Eq. (la), precisely such a solution expansion has been described in detail by Leung, 
Tong and Young [5] . 

The natural correspondence-principle version of the relativistic free-particle Schrodinger equation was 
iterated by Klein, Gordon and Schrodinger for no physically motivated reason, but merely in an effort to rid 
it of its calculationally unpalatable square-root Hamiltonian operator [6, 1, 7]. If this iterated equation is 
still regarded as a complex-valued quantum-mechanical entity, a large class of completely extraneous, highly 
unphysical unbounded-below negative-energy solutions are injected by that iteration. These also destroy its 
probability interpretation, and the fact that it depends on only the square of a Hamiltonian cuts it adrift 
from the Heisenberg picture and Ehrcnfcst theorem. However, if this iterated equation is regarded as the 
description of a classical, real-valued field, it thereupon becomes strongly analogous to the classical wave 
equation, and has an eminently sensible nonnegative conserved energy [3, 1]. This classical Klein-Gordon 
equation is as well one of those classical equation systems which is linearly equivalent to a Schrodinger 
equation: it quite marvelously chooses to be equivalent to precisely the Schrodinger equation with the 
natural correspondence-principle square-root Hamiltonian operator which Klein, Gordon and Schrodinger 
had tried to sideline by concocting it. 

It is a pity that Klein, Gordon and Schrodinger had no idea of the theorem presented by this paper, 
and thus were not equipped to unearth this astonishing fact themselves. If they had but grasped the 
full consequences of the real-valued classical Klein-Gordon equation, they might well have abandoned their 
physically unmotivated rejection of the correspondence-principle mandated relativistic free-particle square- 
root Hamiltonian operator (|cp|^ -|- m^c*) = [7, 1]. 
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